Macroscopic vs. microscopic models of circular Quantum Dots 
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The confinement mechanism of electrons in gated circular quantum dots is studied in a sequence 
of models, from self-consistent 3D Hartree calculations to the semiclassical model of Shikin et al. 
Separation of the vertical from transverse confinement allows accurate 2D Hartree calculations. For 
moderate size dots containing ~ 50 electrons, 2D Thomas-Fermi can also be accurate. Finally, a 
Shikin type parameterization of the electron density allows quantitative study of the importance of 
the many contributions to the confinement potential. 
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I. INTRODUCTION 

Quantum dots defined by electrostatic confinement 
continue to be an important tool for exploring the physics 
of electrons on the mesoscopic scale. Numerous Coulomb 
blockade experiments have provided information on level 
spectra and capacitances. At present, more detailed 
knowledge is becoming available, since new experimental 
methods allow mapping of the charge distribution inside 
such dots and even the details of the individual electron 
wavefunctions can be resolvedi'^'^'^. 

The most accurate electron wave functions are pro- 
vided by 3D Poisson-Schrodinger simulations^!^. These 
methods are highly developed but are computer inten- 
sive and therefore not well suited for systematic stud- 
ies over a wide range of parameter values^. Further- 
more, numerical results are seldom conducive to devel- 
oping intuition about the importance of specific param- 
eters. In parallel with these 3D approaches, a num- 
ber of simpler two dimensional (2D) models have been 
developediS'iiii^'i^ii^'i^ . One of their key assumptions is 
strict confinement of the electron gas to a plane (2DEG), 
and often the transverse confining potential in that plane 
is taken arbitrarily to be a parabola or other simple an- 
alytic shape. Energy density functionals have been used 
to further simplify the computations while still includ- 
ing the main terms in the Hamiltonian. Free parameters 
available in such models can be adjusted to fit the data 
from the experiment being analyzed. The advantage of 
2D models is that they require much less numerical effort, 
simplify the interpretation of results and allow system- 
atic analysis. On the other hand, empirical parameters 
of these models may conceal or distort some aspects of 
the physics at play. 

Our primary objective is to connect these two ap- 
proaches in a detailed and quantitative manner. To do 
so with good accuracy and a reasonable amount of com- 
putation, we restrict our attention to axially symmetric 
dots, which has the advantage of reducing the number of 



degrees of freedom. 

The article is organized as follows: in Section II we 
describe the specific heterostructure layout that will be 
studied. We construct the total potential acting on the 
electrons, including the Coulomb interaction (direct and 
mirror terms), and the electrostatic confining potential 
due to the gates and donors. The resulting Hamiltonian 
is solved numerically in the 3D Hartree approximation. 
By varying the gate voltage we present results for dots 
of different sizes. Our first step towards simpler models 
with a more transparent interpretation is the factoriza- 
tion ansatz between longitudinal and transverse compo- 
nents of the electron wavefunctions. This describes the 
dot in an effective 2D Hartree approximation. We end 
this section by showing that the predictions of this ansatz 
are in very good agreement with the full 3D Hartree cal- 
culations. 

Section HI is devoted to the construction of effective 
confining potentials and coulomb interactions for use in 
2D models, which incorporate finite thickness effects. We 
show in particular that there is a substantial difference 
at short distances between the effective Coulomb interac- 
tion and the l/r law. In Section IV we define an energy 
density functional based on the Thomas-Fermi (T-F) ap- 
proximation for the kinetic energy and, via a variational 
prescription, find equations to determine the correspond- 
ing electron density, along with the coulomb, confinement 
and kinetic energies. 

In Section V, we show how to improve on the semi- 
classical model of Shikin et alr^^^. We do so by taking 
into account the semiclassical kinetic energy, the mir- 
ror coulomb terms and the role of finite thickness of 
the 2DEG on the Coulomb interaction. We replace the 
Thomas-Fermi densities in the variational equations by 
the simple parametrization first introduced by Shikin et 
al}^'^'' and later extended by Ye and Zaremba^^. We 
show that after these simplifications, the numerical effort 
is reduced to a minimum: solving a very simple nonlinear 
equation for the dot radius R. The resulting predictions 
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for bulk properties such as the number of electrons, root 
mean square radius, kinetic, confining and Coulomb en- 
ergies are quantitatively accurate. Finally, by restricting 
the parametrized density to the simplest Shikin form, but 
keeping all relevant potential contributions, wc show that 
the dot size systematics can be obtained from an analytic 
expression where each term has a clear interpretation and 
expresses the physics quantitatively. 



II. MICROSCOPIC MODELS: 
SCHRODINGER. 



POISSON- 



A. The heterostructure 

Fig. [1] shows a section of the heterojunction that we 
will use for reference. It consists of planar layers grown 
on a GaAs substrate which we model as infinitely thick 
and free of acceptors. From bottom to top there are 1) 
an AlGaAs spacer layer of thickness s , 2) an AlGaAs 
donor layer of thickness d, uniformly doped with donor 
concentration pd, 3) a GaAs cap layer, undopcd and of 
thickness c. The cap layer is covered by a metallic gate 
except for a circular opening of radius 5*0 exposing the 
GaAs cap layer. Since our main purpose is to compare 
the predictions of several models, rather than to produce 
an accurate description of a specific experiment, we in- 
troduce several simplifying assumptions that will make 
our calculations less cumbersome, but which should not 
affect the physics to be explored. The first of these is 
the gate geometry just described. Using expressions to 
be given below, we have checked that for rings of a few 
hundred nanometers width, and same inner radius 5*0, 
the potential due to the gates is practically identical to 
that of the simplified gate used here, at the location of 
the dot. Also for simplicity we neglect the differences 
in dielectric constant £r between the different materials, 
and use a common energy-independent effective mass for 
spacer and substrate. 
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FIG. 1: Section of heterostructure layout. Lengths in nm. 
Since in an earlier worki^ we modelled quantum wires 



using similar microsopic models we have adopted the 
same geometrical and physical parameters for the het- 
erostructure: c = 24 nm, d = 36 nm, s = 10 nm, and 
Pd = 6.10"'' nm-3, Er = 13.1, m* /m^ = 0.067. These 
values are taken from a model heterostructure considered 
by Laux et al."^^ in their pioneering work on wires. We 
have discussed how to deal with partially ionized donor 
layers in previous work^^. Here for simplicity we will 
assume the donors are totally ionized. (This is often 
achieved in experiments by illuminating the sample.) Fi- 
nally, we assume that the Fermi level is pinned at the ex- 
posed surface. This again would depend on how the sam- 
ple has been processed, but for modelling purposes it has 
the advantage of being a definite prescription, which has 
been used frequently in a variety of problems^ i^^i^^i^^i^'* . 
It is also a useful reference for other models of surface 
state o^^'^^ '^'^. As indicated in Fig. [U we choose the 
origin of coordinates at the center of the exposed sur- 
face, and the z axis orthogonal to the layers. Vectors 
from that axis, in planes parallel to the layers, are de- 
noted by letters r or s. Followingi^, we construct the 
electrostatic potential acting on the electrons by adding 
the contributions from a) the gate, b) the ionized donor 
layer, and c) electrons in the dot. The potential due to 
a metallic gate, at voltage Vo and with a circular hole of 
radius 5*0, can be derived from expressions given in Sect. 
3.6 of^i; 
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where e is the electron charge ( e < 0.) Using the bound- 
ary condition for Fermi level pinning. 



eVg{s,Q)^eVo Q{s 
one finds, for z > 0, 
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Here E{k'^) is the complete elliptic integral of the second 
kind as defined in2^. To this one must add the elec- 
trostatic potential due to the donor layer, including the 
corresponding mirror term which maintains the bound- 
ary condition at the surface. As ii>i^, in the substrate 
and spacer layers this contribution reduces to a constant 
additive term: 



eVd{z) 



-pd d{c+ -) 



(4) 



The third contribution is due to the Coulomb interaction 
between the electrons. Including the mirror terms it is: 



eVe{f,z;r',z') 
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y/\f-r'\^ + {z-z'y 
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Appropriate band offsets in the AlGaAs layers, eVbo = 
0.23 eV, are included, but for brevity we do not write 
them in the expressions below. Choosing the Fermi level 
as the origin of energies and denoting by eVg the binding 
energy of the surface states with respect to the conduc- 
tion band edge, the latter will be located at: 

eVg+d{r, z) = eV, + eVg(r, z) + eYd{z) . (6) 

We have taken eVg = 0.7 eV. 

The Hamiltonian for N electrons in the dot is: 



i=i ^ 



N 

+ ^ eVein,zf,rj,Zj) , (7) 
i<j=i 

which we solve in the Hartree approximation. 



B. 3D Hartree approximation 

We take a trial product wave function 

N 

^{fi,Zi;....;fN,ZN) = Y\_ipi{ri,Zi) , (8) 

i=l 

and minimize the trial energy 

Eh =< > , (9) 

under constraints that impose orthonormality of the or- 
bitals ipi. To guarantee equilibrium with the surface 
states, all single particle states of negative energy are 
occupied. The variational method gives 



_2m* 
Z7_H-(f, z) 

with 



+ eVg+diri.Zi) + UH{'ri.,Zi) 



ipi = Siipiin, Zi) 



df'dz' eVe{f,z;r',z') pif',z') 



N 



(10) 



(11) 



We will consider only dots with an even number of elec- 
trons. In that case the ground state solution has axial 
symmetry, so that 

p(r, z) = p(r, z) , Unir, z) = Unir, z) , (12) 

and the wavefunctions factorize as 

-u„i(r,z) e*'^ 



V'ni(r, z) 



27r 



(13) 



We have omitted for simplicity the spin component of 
the wavefunction. Furthermore, eq. [10] simplifies to a 
two-dimensional Schrodinger equation 

e^iu^i{r,z) = -—l— + — u^,ir,z) + 



iVg+dir,z) + UH{r,z) + 
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2m* 



Urd ■ (14) 



As usual, we have solved the system formed by eqs. [TUl 
and [HI by numerical iteration. Before presenting results, 
we explore factorized approximate solutions of eq. [TJj . 



C. Factorization ansatz 

It is well known that large confining energies in the 
longitudinal direction support the validity of approxima- 
tions where the z dependence in the wavefunctions is fac- 
tored out from the other degrees of freedom. Guided by 
the success of the ansatz employed for linear wiresi^.^ we 
use it here for dots. We begin by approximating the two 
potential terms in eq. [T3| as 

eVg+d{r,z) + UH{r,z) = 

~ eVg+d{0,z)+m,{r) + UH{r^ (15) 

where the effective two-dimensional confining, AJ/c, and 
Hartree, UH^ potentials are defined as: 



Ai7e(r) 



dzA^{z){eVg+d{r,z) - eVg+d{Q,z)) 



Unir) = / dz A\z)UH{r,z) 



(16) 



with a weight function ^^{z) to be specified below. For- 
mally the range of integration over z extends from — oo 
to -l-oo, but in practice the A{z) are non-zero only in the 
spacer and the substrate. Making these approximations 
in eq. [TOj gives 



P 



2m' 



+ eVg+diO, z) + AC/,(r) + Unir) 



that has solutions 

^^\r,z)^A{z)cj,,{r) 



(17) 



(18) 



with separately normalized A{z) and 0^, which satisfy 
Schrodinger equations with eigenvalues Ez and 



Pz 

2m* 



eVg+diO, z) 



pI 



2m* 



+ AU,{r) + UH{r) 



A{z) 



EzA{z) (19) 
(20) 



Note that we have chosen the previously unspecified 
weight function A{z) to be the common longitudinal com- 
ponent of the wavefunctions in [THj We then have 



E, 
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eP = NE,+Y,<^^\^ + ^U, + -UH \cf>^ (31) 

The numerical process starts by constructing eVg+t;(0, z) 
and solving eq. [191 After that A{z) is fixed and we 
solve the two-dimensional Hartree problem iteratively, 
constructing the potentials from eqs. 1161 and solving 
equation [20] for the (j)i{r), i — 1,...N. For dots with an 
even number of electrons, axial symmetry allows separa- 
tion of radial and angular factors in the 0i(r), as done 
in eq. I13( and eq. [20| reduces to a one-dimensional equa- 
tion for the radial components, which is straightforward 
to solve. 
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angular momentum I 



to filling of the various orbitals. This is displayed in Fig. 
[3| where we have plotted: 

anir) = J p{r,z) dz. (22) 

for N = 40, 54, 72, and the corresponding quantity in 
the factorized approximation, 

N 

1=1 



0.0030 
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FIG. 2: Energies of single electron levels, corresponding to 
Vo = -1.27 eV and So = 210 nm. Circles: 3D Hartree, x 
symbols: factorization ansatz. 



D. Results 

In Fig. [2] we show, for iV = 54 electrons, the orbital 
energies as a function of angular momentum. One sees 
that results using the factorization ansatz (labelled 2-1-1 
D) are in excellent agreement with the exact Hartree 
(labelled 3 D). The various bulk energies are also ac- 
curately predicted by the factorization ansatz: To give 
two examples: 1) when Vq — —1.32 V, we find a 
ground state with = 40 electrons. The kinetic, 
confining. Coulomb and total energies are found to be 
0.720(0.729) eV, -3.263 (-3.272) eV, 1.215(1.214) eV 
and —1.327 (—1.329) eV respectively. The quantities in 
parentheses correspond to the ansatz. 2) when Vq = 
-1.27 V, iV = 54, the respective values are 0.972 (0.985), 
-5.076 (-5.095) , 1.947 (1.949) and -2.160 (-2.161) eV. 

A 3D plot of the electron density p{r, z) of the Hartree 
approximation shows that the dot extends about 20 nm 
along the z axis with a smooth variation very similar to 
the square of the Airy function. This will be exploited in 
Appendix A. Along r the density shows oscillations due 



FIG. 3: 2D electron densities (j{r) computed in the Hartree 
approximation (continuous lines) and the factorization ansatz 
(dashed lines). A*' = 40, 54 and 72 electrons. 

Again, the agreement between the Hartree (3D) and 
the factorized densities is very good. A plot of the same 
densities of Fig. [3| on a logarithmic scale, shows that 
the agreement for the large r is excellent down to val- 
ues as low as anir) — 10^^ nm^^. This suggests that 
the factorization ansatz might be accurate in studies of 
tunneling across a ring of finite width. The root mean 
square radii are also accurately predicted: for = 40 
electrons we find 57.5 (57.6) nm, whereas for iV = 54 the 
two values agree to three digits, 63.5 nm. 

Fig. [H shows the Hartree and confining potentials. 
Both are nearly parabolic for the radii corresponding to 
the area occupied by the dot (?■ < 100 nm), and their 
sum is almost constant in this range of r. We also find 
that the kinetic energy is small on the scale of the figure. 
These findings will be used in Sections IV and V to sup- 
port, and improve on, the macroscopic model of Shikin 
et al^. 



III. POTENTIALS IN THE 2D LIMIT 

Here, our aim is to develop the relation between the 
3D calculations of the last section, and purely 2D models. 
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FIG. 4: 2D potentials A[/c(r) (dashed line), C/j/(r) (dotted 
line) and their sum (continuous line). A'^ = 54, corresponding 
to Vo = -1.27 eV and So = 210 nm. 



We make use of a separable expansion of the confining 
potential developed in Appendix B: 



eVg{Q,z) + h,{z)r^ + ]^ki{zy ... (24) 



eVg{r,z) 



Expressions for the coefficients are given in eq. IB5I To 
have a workable expression for A[/c(r) defined in eq. [TBI 
we approximate the integral over z by replacing the exact 
A{z) by Airy functions (see Appendix A). This leads to 



(25) 



with the functions kdz) and k4^{z) replaced by their val- 
ues ai z = z given in eq. IA12I From now on we will use 
this truncated form of the confining potential, U2D,c{'r) 
with terms parabolic and quartic in r. 

Guided by eq. [121 we define the effective 2D coulomb 
interaction as 



A^{z)A^{z') dz dz' 

y/\f^f'\^ + {z-z'y ' 



(26) 



We now look for a suitable approximation that simpli- 
fies the averaging over the A{z). To do so we again use 
that the A{z) are well approximated by the Airy func- 
tion. With notation introduced in Appendix A, we write 



^ V2Dis) 



Ai'^{u)Ai^{u') du du' 
■y/ s^a^ + {u — u'Y 



(27) 

with s = \r — r'\. In this approximation, the ratio 
V2d{s)/oi depends only on the product sa. The pa- 
rameter a contains all the information on the confining 
potential and on the specific heterostructure. We have 
checked that to a very good approximation the effective 
interactions obtained numerically from eq. [26] for differ- 
ent values of and the gate potential Vq also satisfy 



this simple relation. The continuous line in Fig. [5| corre- 
sponds to the curves (e^/e^) (V2_D(s)/a) computed from 
eq. [26|for a range of values of Sq from 150 to 260 nm and 
of Vb from —1.12 to —1.42 V. On the scale of the figure 
they are all superimposed. The upper dashed line is for 
the pure 1/s potential and shows that the finite thickness 
effect (from averaging over A{z)) is substantial. For the 
simplified 2D models that we will introduce below, it will 
be very useful to have a simple analytic approximation 
for the effective interaction: empirically, we have found 
that the expression 



3 VM^) 



v/s2+t2(a) 



(28) 



with t^{a) = l/(2a2) gives a fairly good fit to the long 
range part of the interaction, and is very simple to use. 
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FIG. 5: Effective 2D Coulomb interactions. Continuous line: 
effective interactions defined in eq. 1261 Upper long dashed 
line: 1/s interaction. Lower short dashed line: parametriza- 
tion of eq. [2HI 

Its physics content is also clear: we have replaced the 
(z — z')2 in the denominator of eq. [2S| by a constant, 
t, which plays the role of an effective thickness of the 
2DEG. In this approximation the effective thickness is 
proportional to 7"^/'^, 7 being the slope (in z) of the 
confining potential. In the following we will show that 
this effective interaction reproduces the direct Coulomb 
potential in the plane of the 2DEG very well. 

A similar average has to be performed to define an 
effective mirror interaction. However, the width of the 
A{z) functions is small compared to the value of z -I- 
z' ~ 2zm = D, with z — Zm the maximum of A{z), and 
therefore it is a good approximation to write 



V2D,mi\r-r\) = - 



A^{z)A^{z') dz dz' 
^\f- + (z + z'Y 

^ (29) 
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IV. THE 2D ENERGY DENSITY 
FUNCTIONAL. 

The Hartree approximation is expected to be valid 
when the number of electrons in the dot is not small. In 
that case, the kinetic energies involved are small com- 
pared to confinement and coulomb'^", and a Thomas- 
Fermi energy density functional becomes a sufficient ap- 
proximation when the main interest is in the bulk proper- 
ties of the system, rather than in details of each individ- 
ual electron wavefunction and energy. Besides requiring 
much less numerical effort than solving eq. HD] for all the 
this simpler model will allow us to obtain several ana- 
lytic approximations that clarify the roles of the physical 
parameters of the system. 

We replace the two-dimensional density constructed 
from the solutions of eq. [20] by its Thomas-Fermi ap- 
proximation, aTpir). By definition the latter is to be 
be determined from a variational condition on a suitably 
simplified expression for the total energy. Starting from 
eq. [21] we write 



Et = NE, + Ek^■ 



Ec^ 



Ecoul 



(30) 



with 



N = j cTTpir) df 

f 

Ekin = -^^^ / '^Tpir) df 



Econ = J ^Uc{r) UTpir) df 

Ecoul = {V2D{\f-f'^) + V2DM\f-^\))<^Tpir) 

. aTp{r') dfdf , (31) 

where we have used the well known expression for the 
kinetic energy Ekin of a 2D Thomas-Fermi gas. Defining 
now: 



E{^i) =Et - ^J.N 



(32) 



with a Lagrange multiplier /i to fix the number of elec- 
trons, the variational condition 5E{^) / Sutp leads to 



^i = E, + — 7raTF(r) + AC/e(r) + UH,Tp{r) (33) 

TO* 



where: 



(y2D{\f-r\) + V2D^A\r-r\)) 



o'Tpi'f') df , 



(34) 



These two equations define the iterative Thomas-Fermi 
approach. The process is started by solving the 
Schrodinger eq. [12] for A{z) to determine E^. Then AUc 



of eq. [25] is computed. The choice previously made for 
the origin of energies fixes fi = 0. Then by iteratively 
solving eqs. [33l and [34l we determine (Ttf (and N). Note 
that we do take the longitudinal energy Ez into account 
in a realistic way, and therefore the 2D Fermi level of the 
dot is not but —E^. If we had taken it to be zero the 
results would be completely different. 
We will not present the numerical results found in this ap- 
proach, because they are very similar to those we will de- 
rive next, using parametrized analytic forms of the den- 
sity. These prove to be more convenient for obtaining 
bulk properties in analytic form. 



V. PARAMETRIZED DENSITIES 

A. Extended Shikin model 

The 2D Thomas-Fermi model is based on a semiclas- 
sical approximation for the kinetic energy. Some time 
ago, Shikin et o/.— noticed that the dominant contri- 
butions should be those of the coulomb and confining 
potentials and proposed a more drastic approximation: 
a model where both the kinetic energy contribution and 
mirror term were neglected in writing the total energy of 
the dot. Further they assumed that a parabolic confining 
potential 



AC/f (r) = hr^ 



(35) 



With these simplifications they obtained an analytic so- 
lution for the charge distribution: 



CTpir) = cr2 \/l 



i?2 



(36) 



with CT2 and R given in terms of the parameters defin- 
ing the potentials. Subsequently, Ye and ZarembaF^ ex- 
tended this model by adding quartic terms in the con- 
fining potential, and found the corresponding analytic 
solution for the electron distribution 



CTqir) 



02 1 



r2 \ ^1'^ 



a4 1 



r 

i?2 



3/2 



e{R~r) . 
(37) 



We extend the above models by including accurate ap- 
proximations for the contributions of the mirror poten- 
tial, the thickness correction to the coulomb direct po- 
tential and the kinetic energy. We will show that one 
can still obtain analytic expressions and that our exten- 
sion leads to results for bulk properties of the dots that 
compare very well with the Hartree calculations of the 
previous sections. 

Let us begin by remembering that in the limit of a 
strictly two-dimensional quantum dot and a 1/r interac- 
tion, the direct coulomb potential corresponding to the 
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parametrized densities above is analytic 



U2D,dir) 



2e, 



-R 



tT2 1 



1-T2 



r 



when r < R, 



r 

'~R^ 
'-IP 



(38) 



and with the constants T2 and T4 equal to zero. In Ap- 
pendix C we derive approximations for the direct poten- 
tial which take into account the finite thickness of the 
dot using eq. [251 These approximations lead to the same 
expression, eq. I38i but with non-zero values for the con- 
stants T2 and T4, which are given in eqs. IC4l and lC8l 
The mirror potential can also be expressed in a simplified 
form, which is still sufficiently accurate, using eq. [221 



U2D.m{r) = 



N 



(39) 



where = + i?^/2. Note that we have replaced the 
I r — 1 2 by its average over r' on a uniformly charged 
disk of radius i?„ that contains the dot charge. The cor- 
responding power series expansion is 



2D.m 



(r) 



ErL 



1 



r 

21? 



3 

8^4 



(40) 



where according to eq. [371 the number of electrons is 



N 



3 



0-2 



rO-4 



(41) 



Finally, for the small contribution of the kinetic energy 
term, we use the truncated expansion 



m* 
0-4 



0-2 



- 8R^) 



3r2 3r4 \ 



With the approximations described above, eq. 
now be written as 



(42) 



can 



2D.. 



U: 



2D.d 



2D.: 



+ — = (43) 



and should be satisfied for all r < R. Explicit forms for 
each of these terms have been written in eqs. [511 [3H1 SOI 
and . Equating the coefficients at each order of the 
power series, we arrive at 
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(44) 



which is a system of equations in three unknowns: (T2 , 
(74 and R. Since the equations are linear, one can elimi- 
nate (72 and (74 leaving a single nonlinear equation in i?, 
which is then solved numerically, by Newton-Raphson. 
We have already seen that for each value of Vq one finds 
a corresponding E^- Therefore, for each Vq, eqs. [44l de- 
termine a set (72, (74 and R, and these give N and other 
quantities such as the mean square radius 
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or the various energies 
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In Figs. [6lto[Tniwe compare the results of the present 
model to our Hartree calculations (with the separation 
ansatz). Fig. [SI shows the number of electrons versus the 
gate potential. The prediction (long dashed fine) is so 
close to the Hartree results that the two lines can only 
be distinguished when eVg is close to 1.1 eV. For compar- 
ison we also show the predictions when the finite thick- 
ness corrections, or the quartic terms, are omitted. As 
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FIG. 6: Electron number obtained from the various models: 
continuous line: factorized Hartree. Long dashed line over- 
lapping with it on most of the range: parametrized density, 
(Tq(r), with parameters determined as described in Section 
V. Short dashed line underneath: same but omitting the non- 
zero thickness corrections in the coulomb terms {t2 = T4 = 0.) 
Upper dotted line: same but omitting the quartic term in the 
confining potential (^4 — 0). 



expected, the quartic terms are relevant only for large 
dots, whereas the finite size correction is sizeable even 
for small ones. Fig. [7] shows the evolution of the various 
energies per particle: Coulomb (direct plus mirror), con- 
fining and kinetic, as eVb is varied. Again the agreement 
is excellent, confirming the validity of the model for bulk 
properties. The agreement is even better for the fields: 
for the example shown in Fig. 21 those given by the 
parametrized densities are so close that they cannot be 
distinguished from those plotted in the figure. This is a 
very useful result, since one can use them as starting val- 
ues for self-consistent solution of the Hartree equations, 
and convergence is then very fast. 

Fig. [8] shows the root mean square radii. Again 
the variation with Vq is fairly well reproduced. But 
one observes a small systematic underestimate that we 
now examine for a specific case. Fig. [TU] compares 
the radial distributions of charge for a gate voltage of 
Vo = —1.22 V. As expected, the quantal oscillations of 
the Hartree density are more pronounced than those from 
the parametrized density. However, the most significant 
disagreement occurs at the surface, where the form of 
only allows an interpolation between the oscillations of 
CTij. This is clearly shown by their difference, the dot- 
ted line in that figure. Still it is seen that most of the 
charge is correctly located, and this is why the energies 
and fields are so accurately predicted by the model. 

Fig. [5] is a more detailed plot of relevant energies v.s. 
Vq. We show separately the contributions from the di- 
rect and mirror Coulomb terms. The mirror term is quite 
large and should not be ignored in any 2D calculation. 
We also show three comparatively small terms: the ki- 



FIG. 7: Energies per electron v.s. gate voltage. From top to 
bottom: Coulomb, confining and kinetic. Continuous lines: 
Hartree with factorization ansatz. Long dashed lines: predic- 
tions of parametrized densities of eq. |37] and expressions [46] 
to 1471 Short dashed lines: analytic approximations: eqs. 1541 

netic energy (positive) , the Coulomb exchange (negative) 
and the correction to Coulomb due to the finite thickness 
of the 2DEG. The latter correction is obtained from the 
expression for the Coulomb energy in eq. 27] as 

^EcouUd = EcouLd - EcouLd{T2 = 0, = 0) . (48) 




FIG. 8: Root mean square radius v.s. gate voltage. Con- 
tinuous line with signs: Hartree, (factorization ansatz). 
Dashed line with circles: parametrized density with two 
terms. Dotted line with diamonds: parametrized density with 
one term, eq. 1511 Dash-dotted line with boxes: prediction of 
eq. m 

For comparison, the Coulomb exchange energy has also 
been computed following Zarembai^, although the ef- 
fect of the corresponding effective potential has not been 
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included in our model. Notice that its contribution is 
comparable to the effect of the finite thickness on the 
Coulomb energy, a term that pure 2D models neglect. 

6 I , , 1 

5 ^x.^ 

4 . ""^-^ 

3 " ^ 








-1 
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1.15 1.2 1.25 1.3 1.35 

eVo ( eV ) 



FIG. 9: Energies relevant to the equilibrium configurations 
v.s. Vq: from top to bottom: Coulomb direct, Coulomb direct 
plus mirror, confining, kinetic. Coulomb finite thickness effect, 
Coulomb exchange and Coulomb mirror terms. 



where Cfe = Ti^ / {2m*)er/ {-Ke^R) ~ 1.65/i? is the smaU 
correction due to the kinetic energy. Similarly the fi- 
nite thickness corrections are only a few percent. Note 
also that in both cases there is a significant cancellation 
between the corrections in the numerator and in the de- 
nominator, so to a first approximation they cancel. By 
neglecting the mirror and the quartic term contributions 
one arrives at Shikin's simple result 

i?| = -|i. (52) 

Due to other cancellations that we will discuss below, Rs 
is fairly close to the value of R found by solving eq. [HI] 
without simplification. This is shown in Fig. [S] To give 
a specific example, when Vo = —1.27 V, the value of R 
from eq. 2?] is 98.0 nm, the value from eq. [51] is 96.5 
nm, and from eq. [221.^5 = 95.8 nm. One can thus find 
a good approximate solution for R by inserting Rs for 
R on the rhs of eq. [HTJ This gives R = 96.6 nm, and 
on the scale of Fig. [5] the exact and the approximate 
solutions to eq. [5T]are indistinguishable. The agreement 
with the Hartree results is also quite good. This clearly 
shows that eqs. [5T] and [SH can be used to quantify the 
effect of the different terms on the size of the dot. 



B. Analytic approximations for the radii and bulk 
energies 

The values of (74 are typically about ten percent of 
those of (72 • Therefore, if one is willing to trade some 
accuracy for simplicity, (74 can be neglected. In this ap- 
proximation, aq{r), eq. [37l is replaced by <Jp{r) of eq. [36] 
and the direct coulomb potential reduces to: 



2D,d 



2er 



i?CT2 1 — r2 



1 T2 \ 

2(3r2-l)-^-y-^l 



(49) 

which, due to the smallness of T2, is practically parabolic. 
It is clear then that imposing eq. [43|for all r < i? will lead 
to unphysical results. What we do instead is to impose it 
aj as a weighted average over all r, with ap(r) as weight 
function and b) at r = 0, where the quartic terms are not 
relevant. We then arrive at two equations: 



= IkcR^ + ^kiR'^ 
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(50) 



which can again be reduced to a single non-linear equa- 
tion for R: 



R-" = -E, 



1 17^ 4 _ 
1 - -T2 - -- 



71? 



56 



(fc, + |fc4i?2) (1 - T2 



+46) 



(51) 




FIG. 10: Radial distribution of charges: Densities weighted 
with a factor 2nr, for Vo ~ —1.22 V. Continuous line: Hartree 
with factorization ansatz. Dashed line: parametrized density, 
eq. 1371 Dotted line: difference. 

Expanding the second denominator in eq. [5T] and ne- 
glecting small terms gives 



R' 



1 



10 
7 



■T2 



AkiRs 

7 kr 



-6 



(53) 



which is less accurate, particularly for large dots. But, 
the usefulness of this simplified expression is that it dis- 
plays very simply the expected role of each contribution: 
the finite thickness correction and the confining quartic 
term decrease the radius, whereas the mirror term in- 
creases it. For our chosen heterostructure, cancellations 
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between the contributions to eq. [53] leave the final dot 
size very close to Rs- The dependence of Rs on Vq can be 
brought out using analytic expressions from the appen- 
dices: the numerator of eq. 1521 Ez, can be approximated 
by eq. IA21 whereas the denominator is linear in Vq as 
shown in eq. IB5I From these expressions, the predicted 
values of Rs fall between the continuous line in Fig. [5] 
describing the Hartree results and the dashed line corre- 
sponding to solutions of eq. [5TJ Summarizing: due to 
mutually compensating effects, the main trends in the 
dependence of the dot size on the gate voltage can be 
easily understood with the help of the expressions given 
in the appendices for the numerator and denominator of 
eq. [52l A similar analysis can be made for the correc- 
tions to (72 using the second of cqs. [50l Then from eqs. 
im and ST] we find 
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.(54) 



Setting R — Rs and using the analytic expression for E^ , 
eq. IA2I the predictions are as shown in Fig. [TJ they are 
quite accurate, so that one can see by inspection of eqs. 
[S3] the role of each term in determining the bulk energies 
given by the 3D Hartree calculations. 



VI. 



SUMMARY AND CONCLUSIONS 



Starting from a 3D Poisson-Schrodinger model, we 
have studied a sequence of approximations. The 3D 
Hartree results show that although the electrons are con- 
fined around the junction plane, they are spread longi- 
tudinally more than 20 nm into the substrate. We have 
introduced a factorization ansate that is very accurate, by 
factorizing out a common longitudinal component of the 
wavefunctions, as a first step towards pure 2D models. 
In the second half of the manuscript we have focussed on 
improving the semi-classical model of Shikin et al^^. We 
have formulated a variational approach where in addi- 
tion to the classical ingredients, (confinement and direct 
Coulomb potentials), we take into account the kinetic 
energy via the Thomas Fermi approximation, the mir- 
ror coulomb terms and the finite thickness of the 2DEG. 
We have shown that using parametrized densities in the 
variational equations leads to good values for the number 



of electrons, mean square radius, kinetic, confining and 
Coulomb energies. Furthermore, the corresponding total 
potential seen by the electrons is so close to that found in 
the Hartree factorization ansatz, that it provides a very 
efficient starting point for full Hartree calculations. In an 
additional simplifying step we have shown that choosing 
the simplest parametrization for the density, but keeping 
all potential contributions, the dependence of the bulk 
properties on the gate voltage can be described analyti- 
cally. These expressions clarify the role of each parame- 
ter in determining the dot properties. This detailed and 
quantitative connection between 3D simulations and the 
simplest 2D models provides an essential tool for dis- 
cussing the systematics of bulk quantum dot properties 
in a transparent and intuitive way. 
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APPENDIX A: APPROXIMATE 
LONGITUDINAL WAVEFUNCTIONS 

1. Airy functions 

Here we describe how we approximate the longitudinal 
wavefunctions A{z), by Airy functions. It is well known 
that the potential term eVg+d{0, z) is close to linear in 
the substrate, and a steep barrier in the spacer layer. 
Therefore, in the substrate we approximate eq. [TH] by 

:A['{z) + 7(z - zj)Mz) = EiMz) (Al) 



2m 



where zj is the location of the junction between substrate 
and spacer, and 7 a constant slope whose value will be 
fixed below. In this approximation: 

47re^ d z 7 

Ez = Ei+ eVs Pdd{c +-) + eVo ^ . 

(A2) 



Introducing the change of variables 

u = a{z - zj) + ^0 , 
with the identifications 



(A3) 



El = -^07'/' 



2m* 



1/3 



(A4) 
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eq. [XT] reduces to the Airy equation in the form given on 
page 446 of Abramowitz and Stegun3^ 



A'l{u) - uAi{u) = 



(A5) 



In the hmit of an infinitely steep spacer barrier, Ai{z) 
must vanish at z = zj so the first zero of the Airy func- 
tion Ai{u) must occur at this point. Since the location 
of that zero is at = —2.33810, the form of Ai{z) is 
completely determined. 



Ai{z) = MM {a(z - zj) + Co) e(z - zj) 



(A6) 



The normalization, and the location of the maximum are 
also analytic 



= zj + -(eo-^o) , 



(A7) 



with = —1.01879. To have a good approximation for 
the potential where the wavefunction is largest, we choose 
to define the slope 7 as 



2. Averages over z 



As will be shown in Appendix B, eq. IB5[ the confining 
potential can be expanded as a power series in . To 
extract from it an effective 2D confining potential using 
eq. [inione has to perform an average over z with weight 
function A^(z). We approximate this average using the 
Ai{z). For the quadratic term one has to determine: 



0C.2 



A^{z) dz , 



(A9) 



and in this expression the denominator varies little 
(z^ << Sq) in the range where A{z) is nonvanishing. 
We remove it from the integral, leaving an average 

zA^(z) dz ~ zj + J{z - zj)Af{z) dz = zj - . 

(AlO) 

Combining this approximation with the second of eqs. 
[7] we find: 



kc - Uz) - ^eVo (^2 +^°2)5/2 ' 



(All) 



7 = 



d{eVg+d{0,z)) 



dz 



02 
^0 



(A8) 



For a given Vq, and the A{z) determined by solving eq. 
1191 we construct its analytic approximation by shifting 
Ai{z) to make the maxima of the two coincide. In Fig. 
miwe compare the shifted Ai{z) to the numerically deter- 
mined A{z) for a representative gate voltage. One sees 
that the Ai forbid tunneling into the spacer, but that 
their shape in the substrate is correct. 



with 



Co ~ "^o _ 
a 3 a 



(A12) 



where Zm is the location of the maximum of A{z). The 
quartic term in the confining potential can also be ap- 
proximated with the same prescription, neglecting the 
variation with z in the denominator of the term in eq. 
IB5I and keeping only the dependence on z in the numer- 
ator. 
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FIG. 11: Longitudinal wavefunction as determined using the 
factorization ansatz (continuous line), compared to the Airy 
approximation described in Appendix A (dashed line.) Vo = 
-1.27 V. 



FIG. 12: Confining potential {Vg{r,z) - Vg{Q,z))/Vo (solid 
line) compared to the power series expansion in with one 
(long dashed line) and two terms (short dashed line ). z = 70 
nm in this example. 
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APPENDIX B: POWER SERIES EXPANSION 
OF THE GATE POTENTIAL 



APPENDIX C: FINITE THICKNESS 
CORRECTIONS 



From eqs. [T] and [21 when z > we have 



eVg{r,z) 



f2TT 

s ds I dO 



27r Jso " ^0 (z2 + |f_ 5^2)3/2 ' 

(Bl) 

To determine the series expansion in powers of r for 
fixed z, we use an auxiUary expansion that generates the 
derivatives d"{eVg)/dr'^ in a straightforward way. We 
expand 
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and integrate over 9 term by term. To illustrate the steps 
we write only the first three terms in the expansion. The 
integration over 9 gives 



eVg(r, z) = eVbz / s ds 
J So 
15 r2s2 
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(^2 + ^2)3/2 2 (z2 + s2)5/2 



15 
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These integrals are analytical, so that 
eVg{r,z) = eVo 
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To get all the contributions of order r^, one has to include 
two more terms in the expansion of eq. IB2I We skip the 
details and quote the full result 



eVg{r,z) = eVo 
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(z2 + 52)l/2 4(^2 + ^2)5/2 
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(B5) 



In Fig. [T^l we show the various results for (Vg (r, z) — 
Vg{0, z))/Vb for z = zj = 70 nm, as a function of r. As 
can be seen, the power series converges nicely towards the 
exact result. We have checked that similar convergence 
holds for other values of z. 



Here we derive the effective Coulomb potential to be 
used when treating the dot as a 2DEG. We first compute 
the exact correction at ?■ = for the two terms of the 
parametrized density, eq. 1371 Then we show that to a 
good approximation a simple parametrization based on 
this correction can be used when r =/= 0. 



1. Potentials at r = 

Using the effective interaction defined in eq. [551 
the first term of the parametrized density, the potential 
that we wish to compute requires the integral 



(CI) 



With the change of variable: x = r'^, the integral takes a 
form that can be easily evaluated: 

m = + ^ (i + ;J) (f - ^r^t^^ I) ■ (C2) 

and when t'^/R'^ «1 



so that defining 
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one can finally write the finite thickness correction as 
I{t) = (1 — T2)I{t — 0). We have checked numerically 
that indeed this approximation is quite good for the 
usual values of R and the value of t given in eq. [551 

For the second term of cTg, eq. [371 the integral of 
interest is 
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Making use of integrals found on page 83 o^, we arrive 
at 
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This result reduces to Iq{t = 0) = 37ri?/8 in the hmit 
t = 0, for which the integration is trivial. Then we define 
(4 = Iq{t)/Iq{0), and write 



= 3^ 



3 /tt R-" -t-" 



f\ t 



(C7) 



If t'^/R^ << 1 can be neglected, the expression for ^4 
simplifies to 



C4 = 1 - T4 = 1 - 



10 t 
3^ R 



(C8) 
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FIG. 13: Direct Coulomb potential for the first term of the 
parametrized densities crq(r), eq. |37l Values of corre- 
sponding to Vo = -1.17 V (top) and Vo = -1.27 V (bottom.) 
Continuous lines: exact numerical calculation using eqs. [26] 
and 1341 Dotted lines: zero thickness eq. IC9I dashed lines: 
including the finite thickness correction: eq. ICIOI 



2. The direct Coulomb potential when r 7^ 



When the 2D interaction is taken as 1 /r the Coulomb 
potential for the density (7p{r) is analytic: 
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In Fig. [13] we compare this curve, dotted lines, to the 
exact numerical calculation using eqs. [25] and [M] As 
expected there is a sizeable discrepancy where the dot 
charge is large, and as the distance from the center of the 
dot increases the thickness effect decreases. The dashed 
lines are our empirical prescription for extending the fi- 
nite thickness correction to non vanishing r: 
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As shown in Fig. 1131 this prescription accurately repro- 
duces the exact calculations. For the small contribution 
of the second term in the parametrized density we use a 
similar prescription: 
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and the results (not shown) are also sufficiently accurate. 
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